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CN ' Abstract 



This paper is devoted to a proof of optimal regularity, near the 
initial state, for weak solutions to the two-phase parabolic obstacle 
problem. The approach used here is general enough to allow us to 
"^ ■ consider the initial data belonging to the class C^'^. 



1 Introduction. 

^ ■ We consider a weak solution of the two-phase parabolic obstacle problem 

00 ; H[u] = f{u) = A+X{«>o} - A X{u<o} in Cw = Biox]0, 1], (1) 

• ! u = ip on Bio, (2) 

^ ■ and u satisfies some boundary conditions on the lateral surface of the cylin- 

der CiQ. Here H[u] = Au — dtu is the heat operator, A^ are non-negative 
constants such that A"*" + A~ > 0, xe is the characteristic function of the set 

^^ \ E, and Biq = {a; : |x| < 10}. Observe that equation ([1]) is understood in 

H ■ distributional sence. We suppose that a given function ip satisfies 

^ e C''' (i?io) . (3) 

We suppose also that sup \u\ ^ M with M ^ 1. 

Cio 
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For a C^ n C°-function u defined in Cio we introduce the following sets: 

n^(u) = {{x,t) e Cio : ±u{x,t) > 0} 
A(u) = {{x,t) e Cio : u{x,t) = \Du{x,t)\ = 0} 
T{u) = d{{x,t) G Cio '■ u{x,t) 7^ 0} n Cio is the free boundary. 

We emphasize that in the two-phase case we do not have the property 
that the gradient vanishes on the free boundary, as it was in the classical 
one-phase case; this causes difficulties. Therefore, we will distinguish the 
following parts of F: 

F°(m) = T{u) n A(m), T*{u) = T{u) \ F°(m). 

From [SUW09J it follows that in some suitable rotated coordinate system 
in R" the set F*(u) can be locally described as Xi = f{x2, ■ ■ ■ ,Xn,t) with 
/ G C^'" for any a < 1. 

1.1 Background and main result. 

In this paper we are interested in uniform L°°-estimates near the initial state 
for the derivatives D^u and dfU of the function u satisfying ([I])-([2D- 

Relative interior estimates were obtained in |SUW09j . The corresponding 
estimates up to the lateral surface were proved in |Ura07] for zero Dirichlet 
data, and in |AU09] for general Dirichlet data satisfying certain structure 
conditions, respectively. Unfortunately the proofs presented in |SUW09] . 
|Ura07j and |AU09j do not work near the initial state. To this end, the 
additional investigation of the behaviour of the solution u close to the initial 
state is required. 

Speaking about regularity up to the initial state, we are only aware of the 
results of [UraOTj, |NPPin] . |Shan8j and |Nys08| . In the papers [NysOSj and 
|Sha08] the authors studied the parabolic obstacle problem near the initial 
state for quasilinear and fully nonlinear equations, respectively. In both 
cases, the estimates of the second derivatives D'^u were not considered, and 
hence, only the gradient Du and the time derivative dtu were estimated. The 
results in |NPP10] are most close to those obtained here. Indeed, the authors 
of [NPPlOj considered the parabolic obstacle problem with more general 
differential operator of Kolmogorov type and established the L°°-estimates 
of D'^u and dtU under the assumption that the initial data (p belongs to the 
space C^'". It remains only to note that the two-phase parabolic problem 
with the initial data ip G C^'" was studied in jUraOTj under the additional 
structure assumption that p vanishes together with its gradient. 

Now we formulate the main result of the paper. 



Theorem 1. Let u be a weak solution of ^-^ with a function ip satisfying 
the assumption ([^. Suppose also that sup \u\ ^ M. 

Cio 

Then there exists a positive constant c completely defined by n, M , A^, 
and (f such that 

ess sup {I-D^mI + \dtu\^ ^ c. 

Ci 

Remark 1.1. The result of TheoremUlis optimal in the sense that we require 
from the initial function ip as much regularity as we want to prove for the 
solution. 

Remark 1.2. The cylinder Ciq is chosen only for simplicity. In fact, the 
problem (I^-^ can be treated in Ci+s for arbitrary 6 > 0. In this case, the 
constant c in TheoremUl will also depend on 6. 

Remark 1.3. From TheoremUl and the well-known interpolation theory it 
follows that Du e Clj^^ (Ci). 

The main strategy used in the present paper foUows. At first we prove 
the estimate of the time derivative. We do this in §2 with the help of regular- 
izations. The next step is to obtain the estimates of the second derivatives. 
The analysis of the second derivatives in §3 is based essentially on the famous 
local monotonicity formula due to L. Caffarelli. 

1.2 Notation. 

Throughout this paper we use the following notation: 

z = (x, t) are points in R*^"^^, where x G M" and t G M^; 

I a; I is the Euclidean norm of x in R"; 

Xe denotes the characteristic function of the set E C R"+^; 

f+ = max{f,0}; f_ = max{— 1>, 0}; 

Br{x^) denotes the open ball in R" with center x" and radius r; 

Br = 5,(0); 

Cr = BrX]0,l]; 

d'Cr is the parabolic boundary, i.e., the topological boundary minus the top 

of the cylinder; 

Qr{z^) = Qr{x^it^) = Br{x^)x]t^ — r'^,t^]. Since our main interest are the 

estimates near the initial state, the radius r in Qr{z^) will be always chosen 

such that t° — r^ = 0. 

Di denotes the differential operator with respect to x^; dt = -^^ 

D = {Di, D2, . . . , Dn) denotes the spatial gradient; 

D'^u = D{Du) denotes the Hessian of u; 



Dy stands for the operator of differentiation along tlie direction v G M", i.e. 

|z/| = 1 and 



DyU = ^ i^iDiU. 



i=l 

II • Wp^E denotes the norm in U'{E), 1 < p ^ oo; 

j ■ ■ ■ stands for the average integral over the set E, i.e., 

E 

J meas {E} 

E E 

^ = ^(|x|) stands for a time- independent cut-off function belonging C^(i?2), 
having support in B2, and satisfying ,^ = 1 in i?i. 

J . It is clear that in the annular B2r{x'^) \ B^{x^) the 

function ^.^ ^^o satisfies the inequalities 

\DirA^)\ ^ c{n)r-\ |Ae.,.o(x)| ^ c{n)r-''. (4) 

For future reference, we introduce the fundamental solution 

f exp(-|x|V4t) , ^ ^ 
— ^ — -7^ — -■, for t > 0, 
G{x,t)=l (47rt)-/2 ' ' (5) 

i 0, for t ^ 

to the heat equation. 

We use letters M, A^, and c (with or without indices) to denote various 
constants. To indicate that, say, c depends on some parameters, we list them 
in the parentheses: c(. . . ). We will write c{ip) to indicate that c is defined 
by the sum ||-DVIIoo,Bio + ||v^||oo,Bio- 

1.3 Useful facts 

For the readers convenience and for future references we recall and explain 
some facts. 

Fact 1.4. Let u he a solution of Equation (|7|j, and let e he a direction in M". 
Then 

H[{D,u)^] ^0 m Cio, (6) 

where the inequalities ^ are understood in the sense of distrihutions. 



Proof. The proof of this assertion can be found in |SUW09j (see also |AU09j ). 

D 

We denote 

I{r,v,z^)= I f \Dv{x,t)\'^G{x-x\f-t)dxdt, 

where r g]0, 1], z° = (x", t°) is a point in R""*"^, a function v is defined in the 
strip M" X [t° — l,t°] , and the heat kernel G{x,t) is defined by ([S]). 

To prove the main theorem, we need the following monotonicity formula 
for pairs of disjointly supported subsolutions of the heat equation. 

Fact 1.5. Let ^ := ^(|a;|) be a standard time-independent cut-off' function 
(see Notation), and let hi, /i2 be nonnegative, sub-caloric and continuous 
functions in C2, satisfying 

hi{0,l) = h2{0,l) =0, hi{x,t) ■ h2{x,t) = m C2. 

Then, for < r ^ 1 the functional 

^{r) = ^{r, hi, h2, e, 0, 1) =^/(r, ^hi, 0, l)/(r, ^/^a, 0, 1) 
satisfies the inequality 

^ir)^^{l) + N{n)\\hi\\lcJMlcr (7) 

Proof. For the proof of this statement we refer the reader to (the proof of) 
Theorem 1.1.4 |CK98j (see also Theorem 12.12 in [CS05] ). D 

Remark 1.6. By rescaling one can easily derive the following modification 
of the local monotonicity formula ^: 

HrXn) ^ HR,(r) + |^I|/^iII|q..(.o)||/^21I1q,,(.o) Wr^R = Vr^. 
Here Cr{x) = ^r,xo{\x\) is a standard cut-off function (see Notation) and 

$(r,Cij) = <l'(r,/ii,/i2,Ci?,^°) = ^I{r,CRhi,z')I{r,CRh2,z^). (8) 



2 Estimate of the time derivative 

For e > we consider the regularized problem 

H[u'] = r{u') in Cg, (9) 

u' = u on dBgx]0,l\, (10) 

u' = ifi, on B,x{0}, (11) 

where /^ is a smooth non- decreasing function such that f^{s) = A"*" as s ^ £ 
and f^is) = — A~ as s ^ e; while cpe is a mollifier of if with the radius 
depending on the distance to dBg such that 

supiv?- v^el ^ e, 

B9 

and M satisfies ([I])- ([2]). 

By the parabolic theory, for each e > 0, the regularized problem fP)- flTT|) 
has a solution u^ with Du^ and 9iM^ belonging to L^ (Cg). 

Lemma 2.1. Lei £ > 0, /ei^ m satisfy (1\)-^, and let u^ he a solution of 

([gj-iny. Then 

sup |n^ — m| ^ e. (12) 

Proof. Setting w = -u^ — m we observe that w\„, = u^^} ~ u, and, conse- 
quently, {w — £:)_,_ |g,^_ = 0. Then Eqs. (I9j) and ([1]) together with integration 
by parts provide for arbitrary t e]0, 1] the following identity 



\f{u)-f\u')\{w-e)^dz= 11 -H[w]iw-e)^dz 



B<)X]0,t] BgX]0,t] 

dx + // iDwl'^dz. 



Ifli^'eh?' 



(13) 





{w>e} 



Taking into account the relations f{u) — f^{u^) ^ on the set {u^ > £:} U 
{u < 0} and {w — e)^ = on the set {m^ ^ e} H {u ^ 0}, we conclude that 
the left-hand side of identity (IT3|) is nonpositive. The latter means that 

supw ^ e. (14) 

Cg 

Replacing in identity ( 1T5]) the term {w — e)_^_ by {w + e)_ and repeating 
the above arguments we end up with 

inf w ^ — e. (15) 

Cg 



Combination inequalities flMj) and flT5|) finishes the proof. D 



We observe also that 

Thanks to condition ([3]) we may conclude that dtu^\ are bounded uniformly 
with respect to e. Moreover, for each small 5 > the functions u^ are smooth 
in the closure of the cylinder Cg.^. 

With the estimate (IT2l) at hands it is easy to check that 

ll^tM^hA ^ c, 

and the latter inequality is also uniform with respect to e. Hence we can 
easily deduce the following result. 

Lemma 2.2. For each small 5 > the uniform estimate 

sup \dtu^\ ^ iVi(n, M, A^, (^, 5) (16) 

holds true for solutions u'^ of the regularized problem /f^)-/ f77]) . 

Proof. We set v = dtu^ . It is easy to see that v± are subcaloric in Cg. 

Now we may apply the well-known parabolic estimate (see, for example. 
Lemma 3.1 and Remark 6 from [NtJllj) and get 



supf-t ^ A^W -r v'j^{z)dz + sup f±. 



which implies the desired inequality flT6|) . D 



Remark 2.3. By virtue of Lemma \2.1[ solutions u'^ of the regularized problem 
/fP))-/ l77]) converge to u as e ^ Q uniformly in Cg. 

Remark 2.4. It follows from Lemma \2.S\ and Remark \2.3\ that the estimate 

snp \dtu\^Ni{n,M,X^, if, 1) (17) 

Cr 

holds true for a function u satisfying (CP-IB'- Here Ni is the same constant 
as in Lemma \2.2[ 



3 Estimates of the second derivatives 

Lemma 3.1. Let the assumptions of TheoremUl hold, let z^ = {x^,t^) be a 
point in Ci, and let R = \n^. 

Then there exists a positive constant N2 completely defined by the values 
of n, M, X^, and the norms of ip such that 

sup \Du{-,t)-Dip{-)\^N2R for te]0,t°]. (18) 

Remark 3.2. It is evident that Qqr{z^) C Cj. 

Proof. First, we observe that the assumption ([3]) imphes Aip G L°° (-Bio). 
Therefore, for almost all t e]0, 1] the difference u{-,t) — ip{-) can be regarded 
in Bj as a solution of an elliptic equation 

A (m(x, t) - ip{x)) = F{x) = /(m(x, t)) + dtu{x, t) - Aip{x). (19) 

Due to estimate f lT7|) . the function |-F| is bounded by the known constant up 
to the bottom of the cylinder C7. 

Thus, for a test function 77 G Wq' {Bqr{x'^)) we have the integral identity 

B6r{x0) B6r{x0) 

We set in ( 120|) rj = {u — ip)Q^^ where C'iR{^) '■= ^3R,xo{\x\) is a standard 
time- independent cut-off function (see Notation). Then, after integrating by 
parts and subsequent application of Young's inequality, identity ( 120|) takes 
the form 

/ \D{u - ip)\'^Cmdx ^ ci / (m - ^Xindx 

BeRix°) BeRixO) 

+ C2 / {u-y^)^\DC3ii\^dx. 
B6r{x0) 

It is easy to see that inequality (IT7|) implies the estimate 

sup 1^ — V^l ^ cR^. 

BeR(xO)x]0,tO] 

Putting together fl2T]) and the above inequality we arrive at 

f \D{u{x,t)-ip{x))\^dx^N2R''^^, te[0,t% (22) 

B3r{x0) 
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Now we observe that for almost all t g]0, 1] and for any direction e G M" 
the difference DeU—Df,ip may be considered as a weak solution of the equation 

A {DMx, t) - D,^{x)) = -D,F{x) 

in By. The well known results (see |LU68j . |GT01j ) applied to the difference 
DeU — Deif yield the inequality 



sup \DeU — Dgifil ^ c 



\ 



\DeU — De'fl'^dx + cR sup |F|. 

BsRixO) 
BsRixO) 



Combining the last inequality with the estimate (122|) . we get (TT8l) and finish 
the proof. D 



Lemma 3.3. Let the assumptions of Lemma \3. 1\ hold. Then 

f [ \D\{x,t)\^CUx)G{x-x°,t°-t)dxdt^N3{n,M,\^,ip)R^, (23) 

Jo JM" 

where Cr{x) := C,r,xo{\x\) is a standard time-independent cut- ojf function (see 
Notation) and the heat kernel G is defined by the formula ^. 

Proof. Suppose that e is an arbitrary direction in R" if Du{z^) = and 
e _L z/, where u = Du{z^)/\Du{z^)\, otherwise. From ([3]), fITS]) and our choice 
of e it follows that 

sup \De(p\ ^ cR. (24) 

B2r{x°) 



According to Fact 11.41 the functions v = {Deu)_^_ are sub-caloric in C2, 
i.e., H[v] ^ in the sense of distributions. Since \Dv\'^ + vH[v] = ^Hlv"^] we 
have 



R2 



\Dv{x, t)\%R{x)G{x - x", r - t)dxdt 

B^RixO) 

1 
2 

B2r{x0) 



(25) 
Iff H[v^{x,t)]CUx)G{x-x°,t° -t)dxdt. 



After successive integration the right-hand side of (125]) by parts we get 
I I \Dv\\lGdxdt^- I (y'^'rG) 



B2r{x0) B2r{x0) 

i?2 



dx 





+ f f y Cfi [dtG + AG] dxdt 



B2i?(:!;0) 

/?2 



B2r(x0) 

=: /i + /a + /a- 
It is evident that due to ^ we have 

\DMxM^ 



+ 11 v^[2CRDCnDG + G\DU\^ + GUACR\dxdt 



h^ I l^^^(^^^0)l cUx)G{x-x',f)dx 



B2fl(^'0) 



where the last inequahty provided by 
Taking into account the relation 

dtG + AG = dtG{x-x°,t^-t) + AG{x-x°,t°-t) = for t < t°, 

we conclude that I2 = 0. 

Finally, we observe that the integral in I3 is really taken over the set 
E =]0,i?^] X {i?2/?(a^°) \ -Bij(x°)}. Therefore, in E we have the following 
estimates for functions involved into I3 



\Gix~x',t'-t)\^c^^^^^^cR-; 
\DG{x-x'',f-t)DC{x)\ <c|G(x-x°,t°-t)|- '^"^ ' 



' R{R^ - 1) 



p UR^-t) 



\R^-tY^^I^ 
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Consequently, 

I3 ^ CR-''-^ II v'dxdt ^ Csnpv' ^ CR\ 



E 
E 



where the last inequality follows from flT8|) and (|23 
Thus, collecting all inequalities we get 



a 



I I \Dv\^ClGdz^CR^ (26) 

B2r{xO) 

Inequalities fl26|) mean that we obtained the desired integral estimate for 
all the derivatives D [Dgu) with e _L z/. Similar estimate for the derivative 
D^ (D^u) follows from (^ and Eq. ([1]). D 

Corollary 3.4. Let the assumptions of Lemma IJ.il hold, and let e be an 
arbitrary direction in M". Then 

$e(i?) = $(i?, iDeu)+, (Deu),, Cr, ;2°) ^ N,{n, M, A±, yp), 

where the functional $ is defined by the formula ^, while Cr{x) := ^_R,a;o(|3;|) 
is a standard time-independent cut-off function (see Notation). 

Proof. The desired inequality follows immediately from the definition ([8]) 
combined with (ESI). □ 



Lemma 3.5. Let the assumptions of Lemma \3.1\ hold. letv = Du{z^)/\Du{z^)\, 
and let e be an arbitrary direction in M" if \Du{z^)\ = and e _L z/ otherwise. 
Then 

IIPen)±|||Q,^(,o) ^ iV5(n,M,A±,^)i?"+^ (27) 

Proof. It is evident that inequalities fITS]) and i^^ imply the estimate flTTI) . 

D 



Proof of TheoremUi Let z^ = {x^,t^) be an arbitrary point from Ci such 
that |u(z°)| > 0, and let e G M" be the same direction as in Lemma [3.51 



Since Deu{z^) = 0, it follows that 

C(n)|L'(Deu)(z°)|^^lim<l>e(r), 



11 



where ^e{r) = $(r, {Deu) + , {Df,u)^,QR,z^) with the functional $ defined by 
the formula ([B]), and Cr{.^) = ^r,xo{\x\)- On the other hand, according to 
Fact 11.51 and Remark 11.61 after that, we have for < r ^ i? := vt° the 
inequality 



Application of Corollary 13.41 and Lemma [3.51 enable us to estimate the right- 
hand side of the last inequality by the known constant. This means that we 
proved in the cylinder Ci the L°°-estimate for all the derivatives D {Dgu) [z^) 
with e ^- V. It is clear that the corresponding estimate of D^ {DpU) {z^) can 
be now deduced from Eq. ([1]). 

So, we establish the L°°-estimates for D'^u{z^) for all points 2° G Q"^ fl 
Ci, and these estimates do not depend on the distance of z'^ from the free 
boundary T{u). Since \D'^u\ = almost everywhere on A{u) and the 
(n + l)-dimensional Lebesgue measure of the set T*{u) equals zero, we get 
the uniform estimate of the Lipschitz constant for Du{-,t) for each t G [0, 1]. 

It remains only to observe that the uniform L°°-estimate of dtu were 
established in fITTj) . This finishes the proof. D 
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